INVARIANT INTEGRABILITY CRITERION FOR 
THE EQUATIONS OF HYDRODYNAMICAL TYPE. 
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Cf"^ ' Abstract. Invariant integrability criterion for the equations of hydrodynamical 

Qs , type is found. This criterion is written in the form of vanishing for some tensor 

Qs 1 which is derived from the velocities matrix of hydrodynamical equations. 



oo 



1. Introduction. 



j=i 



Systems of quasilinear partial differential equations of the first order arise in 
different models describing the motion of continuous media. Special subclass of such 
systems is known as a systems of equations of hydrodynamical type. In spatially 
^ i one-dimensional case they are written as follows 

m 

o 

O . (1.1) ul = y^A'4u)ui, where z = 1, 

O 

Q>^ ' Among the systems (1.1) man can consider special subclass of systems possessing 

the Riemann invariants. These are the systems which can be transformed to the 
diagonal form 

(1.2) ul — Xi{u)ul., where i ~ 1, . . . ,n 

o' 

Y\ ' by means of so called point transformations 

KJf - (1.3) M* = u*(u\ . . . ,u"), where i = 1, 

^'. . . 

C^ ' System of equations (1.1) is called the hydrodynamically integrable system if it has 

the continuous set of hydrodynamical symmetries (or hydrodynamical conservation 
laws) parameterized by n arbitrary functions of one variable. For diagonal systems 
(1.2) with mutually distinct characteristic velocities (A^ ^ \j) one has the well- 
developed theory of integration (see reviews [1] and [2] ) . As it was shown in [2] the 
diagonal system (1.2) is integrable if and only if the following condition is satisfied 
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where di — d/du^ and dj — d/duK Such systems are cahed semi-Hamiltonian 
systems, for the property (1.4) itself in Russian papers the term semihamiltonity is 
used. When the diagonal system (1.2) possess this property it can be integrated by 
means of "generalized hodograph method" (details see in [2]) 

In section 2 of this paper we consider the problem of hydrodynamical integra- 
bility for the systems of equations (1.1) with the velocities matrix AUu) of general 
position i.e. eigenvalues of which are mutually distinct. There we managed to prove 
the following fact: each hydrodynamically integrable system (1.1) with the matrix 
of general position is necessarily diagonalizablc^. Summing up this fact with the 
results of [2] we may state the following theorem 

Theorem 1. System of equations (1.1) with the m,atrix of general position is hy- 
drodynam,ically integrable if and only if it is diagonalizable and semi-Hamiltonian. 

Theorem 1 shows that the study of diagonal equations (1.2) is very important. 
But it doesn't exclude the necessity of study of general equations (1.1). Indeed 
the intcgrability test for the equations (1.1) of general position according to the 
theorem 1 should include 3 steps 

(1) test of diagonaHzability, 

(2) diagonalization by means of the point transformation (1.3), 

(3) test of semihamiltonity (1.4). 

First of these steps is implemented by means of invariant geometrical criterion 
from [3]. This criterion consists in vanishing of Haantjes's tensor derived from the 
velocities matrix of the equations in question. This test was first applied to the 
equations (1.1) in [4]. 

Next step is an actual diagonalization. To pass this step one should calculate 
the eigenvalues of the matrix A*(u), find its eigenvectors, properly normalize them 
and then one should solve some system of ordinary differential equations defining 
the transformation (1.2). Because of this step the total integrability test may be 
absolutely inefficient since the case when the system of differential equations is 
explicitly solvable is very rare event. However if we find such solution the third 
step may have only the difficulties in calculations. 

The presence of nonefficient step in the above integrability test of the equations 
(1.1) is due to the absence of the of invariant criterion for testing the semihamil- 
tonity for these equations. The main goal of this paper is to eliminate this essential 
fault of the theory of such equations. As it was noted in [1]: It was Riemann who 
first recognized that the theory of the equations (1.1) is the theory of tensors since 
the components of matrix A* (u) are transformed as the components of tensor under 
the point transformations (1.3). Therefore it is natural to expect that the semi- 
hamiltonity relationships (1.4) can be rewritten in an invariant tensorial form. In 
the section 4 of this paper we construct the tensor, vanishing of which is equivalent 
to (1.4). So invariant integrability criterion for the equations of hydrodynamical 
type is obtained. New integrability test now is absolutely efficient. It includes two 
steps 

(1) test of vanishing the tensor of Haantjes, 

(2) test of vanishing the semihamiltonity tensor. 

In order to construct the semihamiltonity tensor wc use the theory developed 
by Froelicher and Nijenhuis in [5] and [6]. This theory in brief is given in section 



-"^May be this is not new fact but we couldn't find it anywfiere. 
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3. According to the theory of Froelicher and Nijenhuis each smooth manifold is 
equipped with some Lie superalgebra of tensor fields of type (l,p). Note that 
the paper [6] was written in 1956 but unfortunately it wasn't known to specialists 
since for example in [7] the paper of Martin of 1959 is quoted as the first paper in 
supermathematics. 

Authors are grateful to V.E.Adlcr and I.Yu.Cherdantzcv for the fruitful discus- 
sions. They are also grateful to B.I.Sulcymanov for supporting interest to this 
paper. 

2. Hydro DYNAMICAL integrability. 

Let's take the system of equations of hydrodynamical type (1.1) and let's add to 
it another such system with the dynamics by the variable r 



(2.1) ul ^'^Bj{u)ui, where i = 1, ... ,n 

i=i 

Definition 1. System of equations (2.1) is called the hydrodynamical symmetry 
for the equations (1.1) if the equations (1.1) and (2.1) are compatible. 

Now we shall study the question about the existence and the number of hydro- 
dynamical symmetries for the system of equations (1.1). Let the equations (1.1) 
and (2.1) be compatible. Their compatibility conditions are written in form of the 
relationships 

n n 

(2.2) Y.'^B^-T.BiA^ 

n 

J2{dsA^B-; + d^Apt + d,B;A': + d.BtA^) = 

(2-3) 

- Y.idsB'^A^ + dsB^At + d^A^B^, + d.AlBl) 

s=l 

The relationship (2.2) means that the matrices the systems (1.1) and (2.1) are 
commuting 

(2.4) AB = BA 

Second relationship (2.3) also can be written in an invariant form. In order to do 
it let's contract (2.3) with X^X^ , where X^, . . . ,X" are the components of some 
arbitrary vector field X. 

(2.5) [AX,BX] - A[X,BX] -B[AX,X] =0 
This result can be stated as a theorem. 
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Theorem 2. System of the equations (2.1) is a hydrodynamical symmetry for the 
system (1.1) if and only if for any vector field X the relationships (2.4) and (2.5) 
hold. 

Let the operator field A = A(u) from (1.1) have n mutually distinct eigenvalues 
\ = K (u) . Through Xi , . . . , X„ we denote the frame formed by eigenvectors of 
operator A. The choice of eigenvectors is not unique, there is the gauge arbitrariness 
in n scalar factors 

(2.6) X,(u) -^ /,(u)X,(u), where /, ^ 

For the sake of brevity we introduce the following notations for the Lie derivatives 
along the vector fields Xi, . . . , X„ 

(2.7) L, = Lx, 

Mutual commutators of the vector fields Xi , . . . , X„ are convenient to be expanded 
in the frame formed by these fields 

n 

(2.8) i,X, = [X„X,]=^4Xfc 

fc=i 

Parameters cf — c^ (u) in (2.8) are to be called the structural scalars of the frame 
Xi,...,X„. The term structural constants doesn't suit since c*l depend on the 
point u. 

From the algebra we know that the matrix B is commuting with the matrix A 
having mutually distinct eigenvalues then these two matrices are simultaneously 
diagonalized in the frame Xi, . . . , X„. Therefore any operator B satisfying (2.4) is 
completely defined by its eigenvalues /ii = /ii(u). Note that from (2.5) we have the 
relationship 

[AX, BY] + [AY, BX] - A[X, BY]- 
-A[Y,BX] -B[AX,Y]-B[AY,X] =0 

which holds for two arbitrary vector fields X and Y. Let's substitute X = X^ and 
Y = Xj into the above relationship. As a result we obtain that it is equivalent to 
the pair of sets of relationships. First set is algebraic with respect to the eigenvalues 
lii of the matrix B 

c)ki>^j - \k)lJ.t + c]k{\k - \t)l^j + 

+ c)^{\i - \;j)iJ.k = 0, for i ^ j, j y^ k,k^i 

Second set contains the partial differential equations with respect to Hi 
(2.10) L.fXj ^ Kj ^' " ^' , ioii^j 

Here and everywhere below we use the notations Ay — LiXj in terms of Lie deriva- 
tives (2.7). 
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System of differential equations (2.10) is ovcrdctcrniincd. When it is compatible 
the maximal degree of arbitrariness for its solutions is n functions of one variables. 
Let this degree of arbitrariness be actually realized. Then 

(2.11) A*» =Ai»(/i, •••,/«,") 

Let's substitute (2.11) into (2.9). This leads to the functional dependence for the 
parameters /i(zi), . . . , /„(2:„) from (2.11), which contradicts their arbitrariness. 
Therefore the relationships (2.9) should be trivial. For the structural scalars (2.8) 
this gives 

(2.12) 4--0, ioTt^j,j^k,kj^i 

The relationships (2.12) have the important consequences which arc due to the 
following lemma. 

Lemma 1. The linear operator of the general position A = A(u) is diagonalizable 
by means of the transformation (1.3) if and only if the relationships (2.12) hold for 
the frame of its eigenvectors. 

We give the sketch of proof of this lemma. Operator A(u) is diagonal in the 
frame of its eigenvectors. For A(u) to be diagonalizable by the transformation 
(1.3) this frame should be the coordinate frame i.e. structural scalars of this frame 
should be identically zero cf — 0. The relationship (2.12) provides vanishing most 
of these scalars. One can reach vanishing the rest of these scalars by use of the 
gauge arbitrariness (2.6). 

Because of lemma 1 the further analysis of the compatibility conditions for the 
equations (2.10) becomes unnecessary. For the diagonal systems (1.2) such analysis 
was done by S.P.Tsarev in [2]. Note only that as result of such analysis wc add to 
(2.12) the following relationships 



r. (^izi^Vi/ ^'' 



(2.13) 



cL-Aiii. 

= 



c^jAjfe cljKk 



Aj ^ Afe Ai — Xk 



which hold for i y^ j,j y^ k,k j^ i. The relationships (2.13) are the same as 
the semihamiltonity relationships (1.4) but written in the frame of eigenvectors of 
diagonalizable operator A(u). The above considerations prove the theorem 1 in 
the following form. 

Theorem 3. System of the equations (1.1) with the matrix of general position 
possess the continuous set of hydrodynamical symmetries with functional arbitrari- 
ness given by n functions of one variable if and only if it is diagonalizable and 
semi-Hamiltonian. 

Now let's study the similar question about the conservation laws for (1.1). On 
the set of their solutions we define the integral functionals of the following form 

(2.14) F^ f f{u)dx 
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Functional (2.14) is called the hydrodynamical conservation law or the first integral 
for the equations (l.l)ifi^ = when time derivative F is calculated according to 
the dynamics given by (l-l)- This derivative is the following integral functional 



(2.15) G = F = 



P / n n \ 



From vanishing the functional (2.15) we get the vanishing of its variational deriva- 
tives 

J=l fe=l 

This leads to the following relationship for the density /(u) of the primary func- 
tional (2.14) 

n 

(2.16) ^(a.(9fe/4) - d.idkA'i)) = 

k=l 

It is equivalent to the existence of the function T(u) such that 

(2.17) G =fJ2T.{9^fA)ui)dx = / gdx 

Because of (2.17) the condition (2.16) is exactly the condition of vanishing the 
functional G = F = 0. In order to write the relationships (2.16) in an invariant 
form let's choose two arbitrary vector fields X and Y. After contracting (2.16) with 
X^ and Y^ we can write the result of such contraction through the Lie derivatives 

(2.18) (ix-^AY - -^Y^AX - iA[X,Y])./ = 

This result can be stated in form of the following theorem. 

Theorem 4. Integral functional (2.14) is hydrodynamical conservation law for the 
system of equations (1.1) if and only if for any choice of vector fields X and Y the 
equations (2.18) hold. 

The relationship (2.18) is the system of differential equations with respect to 
the unknown function /(u). One should investigate it for the compatibility. Let's 
denote Lif — ipi. Substituting frame vectors Xi, . . . , X„ for X and Y into the 
relationships (2.18) we get the following equations 

n 

(2.19) L,^, -^^Sf^^fc, for^^j 

fc=i 

where the functions Bf^ are defined by the formulae 

(2.20) B^. ^ 4^^-^ ^ ^^ - ^^ 



Xj — Xi Xj — Xi Xj — Xi 
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The equations (2.19) are analogous to the equations (2.10). When they are com- 
patible their solutions have the arbitrariness in n functions of one variable. Let 
such arbitrariness be actually realized. We look for the differential consequences of 
the equations (2.19). Among them we find the following relationships 

n 

(2.21) BJ,L,(^, - BlL.ip, - c%U^k = - ^ Rl,^, 

which hold for i ^ j^j ^ k,k ^ i. The values of i?^- for i ^ j,j ^ k,k ^ i are 
calculated according to the formula 



(2.22) 



n n 

s^i s^k 



The derivatives hiipi, hjipj and hk^k aren't defined by the equations (2.19). This 
gives the arbitrariness in n functions for the solutions of (2.19). When they are 
nontrivial the relationships (2.21) define the functional dependence between these 
derivatives. Therefore they diminish the degree of arbitrariness. In case of maximal 
arbitrariness the relationships (2.12) should be trivial 

(2.23) i?i, - B^, - c^- = 0, ioTt^j,j^k,k^i 

From (2.23) due to the lemma 1 we get the diagonalizability of the operator A(u) 
by means of point transformations from (1.3). Due to (2.20) the equations (2.19) 
are rewritten as follows 



_ ^ji^i — ^ijVj , J 



(2.24) U^, = - -^-'-^-^ + 4^„ for ^ + j 

The compatibility conditions for (2.24) are defined by the quantities from (2.22) 
as i?^- = 0. On taking into account (2.23) these compatibility conditions are 
exactly coincide with (2.13). In spite of the fact that the equations (2.10) and 
(2.24) are different their compatibility conditions are the same and have the form 
of semihamiltonity condition written in the frame of eigenvectors of the operator 
A(u). The above considerations prove the following version of the theorem 1. 

Theorem 5. System of the equations (1.1) with the matrix of general position 
possess the continuous set of conservation laws parameterized by n functions of one 
variable if and only if it is diagonalizable and semi-Hamiltonian. 

3. The Froelicher-Nijenhuis bracket and the Lie 
superalgebra of vector- valued differential forms. 

Let A be the tensor field of the type (l,p) and let it be skew symmetric in co- 
variant components. Then A defines the vector- valued p-form A ~ A(Xi, . . . , Xp). 
Here Xi, . . . , Xp are some arbitrary vector fields. Let B be the second tensor field 
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of the type {l,q) which define vector valued g-form B(Xi, . . . , X^). We shall call 
A and B the vector fields of rank 1 if they are of the following form 

(3.1) A = a«)a B = b(X)/3 

where a and b are vector fields while a and (3 are the differential forms. For 
the tensor fields of the form (3.1) we define the pairing {A, B} (it is known as 
Froelicher-Nijenhuis bracket) 

{A, B} = [a, b] (g) a A /3 - a L^a A /3 + b (g) a A £a/3+ 

^ ■ ' +{-iyaiE>LhaAdf3+{~l)Ph^daAL^f] 

Via ta and tb in formula (3.2) we denote the differentiations of substitution. For 
the r-form uj and for the vector field c the expression LcI-^ is a r — 1-form 

LcUjiXi,. .. ,X,._i) = rw(c, Xi,.. .,Xr_i) 

The operation tc is also known as inner product with respect to the vector field c 
(see [8]). 

Theorem 6. The bracket {A,B} defined for the tensor fields A and B of rank 1 
by the formula (3.2) is uniquely continued for the arbitrary tensor fields of the types 
(l,p) and (1,9) skew symmetric in their covariant components. 

PROOF. Each tensor field of the type (l,p) can be written as a sum of tensor 
fields of rank one as follows 

(3.3) A = ^ A, = ^ ai ai 



The analogous formula can be written for the field B. Therefore the bracket (3.2) 
for the arbitrary A and B can be redefined as 

(3.4) {A,B} = ^^{A„B,} 

i J 

However the expansion (3.3) is not unique. Therefore the definition {A,B} by 
means of (3.4) should be tested for the correctness. The arbitrariness in the expan- 
sion (3.3) for A is defined by the following identities in tensor algebra 

(3.5) / ~^ 

a(8)(Q; + a) = a®a + a(g)Q; 

(3.6) (/a)®a = a® (/a) 

where / is an arbitrary scalar field. The arbitrariness due to (3.5) does not infiuence 
to the value of bracket (3.4) since the relationship (3.2) is additive with respect to 
a and a. Let's ensure that the arbitrariness due to (3.6) also doesn't make the 
influence to the value of {A, B}. In order to do it we calculate this bracket by (3.2) 
first for A — (/a) (g) a then for A = a (g (fa) and after all we compare the results. 
All these calculations are based on the following formulae from [8] 

[fa, b] = /[a, b] - aib/ i/a/3 = /ia/3 + df A l^i3 

Lbifa) == fLi,a + L^fa if^fi = /ta/3 

Since these calculations are standard we did not write them here. Theorem is 
proved D 

Note that forp = g = the bracket (3.2) coincides with the ordinary commutator 
of vector fields. For the arbitrary values of p and q the algebraic properties of this 
bracket are given by the following theorem. 
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Theorem 7. The bracket {A,B} for the tensor fields of rank 1 defined by (3.2) 
and then generalized by (3.4) satisfies the relationships 

{A,B} + (-l)PHB,A} = 

{{A, B}, C}(-l)'-P + {{B, C}, A}(-1)P'? + {{C, A}, B}(-1)^'- = 

because of which it defines the structure of graded Lie superalgebra in tensor fields 
of type (l,m) skew symmetric in covariant components. 

Let A and B be tensor fields of tlie type (f .1) i.e. operator fields. Tensor field 
S = 2{A, B} is tlie vector valued 2-forin. Its values may be calculated by the 
following formula 

S(X,Y) = [AX, BY] + [BX, AY] + 
(3.7) +AB[X,Y]+BA[X,Y] -A[X,BY]- 

- A[BX,Y] -B[X,AY] -B[AX,Y] 

Tensor S is known as the torsion of Nijenhuis for the operator fields A and B (see 

[5], [6] and [8]). 

4. The construction of semihamiltonity tensor. 

First let's recall the classical invariant criterion of diagonalizability for the op- 
erator field A. We mentioned this criterion in section 1 (sec also [3], [4] and [6]). 
Let's consider the particular form of tensor (3.7) 

N = {A,A} 

It is usually called the tensor of Nijenhuis. From (3.7) we obtain 

N(X,Y) = fAX,AYl + A2[X,Y1- 

(4.1) V . y L . J L > J 

^ ^ - A[X,AY] - A[AX,Y] 

Tensor of Haantjes is defined via the tensor of Nijenhuis (4.1) according to the 
formula 

H(X, Y) == N(AX, AY) + A2n(X, Y)- 

(4.2) V , ; \ , J \ , J 

^ ' -AN(X,AY)- AN(AX,Y) 

It has the same type as the tensor N. It is also the vector- valued skew symmetric 
2-forin. 

Theorem 8 (criterion of diagonalizability). The operator A{u) of general position 
with mutually distinct eigenvalues is diagonalizable by means of point transforma- 
tions (1.3) if and only if its tensor of Haantjes (4.2) is identically zero. 

The criterion of diagonalizability in form of this theorem was first proved in [3] . 
It was applied to the systems of equations (1.1) in [4]. 
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PROOF. Because of skew symmetry of bilinear form (4.2) it's enough to test 
vanishing this form only for vector fields X = X^ and Y — Xj from the frame of 
eigenvectors of the operator A(u) with i ^ j- By direct calculations we obtain 



(4.3) H(X„X,)=.^(A,-Afc)2(A, 



2„k 
fc=l 



Because of (4.3) the equality H(Xi,Xj) = is equivalent to (2.12). Then we are 
only to apply the lemma 1. Criterion is proved. D 

Practical use of this criterion for the testing the diagonalizability is based on the 
following formulae for the components of tensors N and H expressing them through 
the components of matrix A(u) 

n 

(4.4) iV^. - Y. i'^tdsA'; - A^d.A^ + A^d.At - A%A^^) 

n n 



(4.5) 



Ag N^j A j AgNj^^Aj + Ng^A^ Aj j 



Let B be the operator field i.e the tensor field of type (1, 1) and let Q be the 
skew symmetric tensor field of type (1,2). Through K we denote the Froelicher- 
Nijenhuis bracket of these two fields K = 3{Q, B}. For the 3-form K we have 

K(X, Y, Z) = B[X, Q(Y, Z)] - [BX, Q(Y, Z)] + 

(4.6) + BQ(X, [Y, Z]) + Q(X, B[Y, Z]) - Q(X, [BY, Z])- 
^Q(X,[Y,BZ]) + ... 

Dots in (4.6) denote 12 summand that can be obtained from explicitly written 
summands in (4.6) by means of the cyclic transposition of X, Y and Z. 

Starting with deriving the invariant semihamiltonity criterion we should note 
that everywhere above (see theorems 1, 3 and 5) the semihamiltonity comes together 
with diagonalizability. It doesn't play the separate role. Therefore we will use the 
following scheme of action: we will construct the tensor vanishing of which gives the 
equations (1.4) after bringing this tensor to the coordinates where the matrix A(u) 
is diagonal. The equations (1.4) are rational. Let's rewrite them in polynomial 
form. In order to do it we introduce the following quantities 

Q^fcij — ^{^i ^ ^j){^i ^ ^k){^j — ^k)dijXk — 

(4.7) - (A, - Aj-)(A, + Aj- - 2Xk)d,XkdjXk+ 

+ (Xi — Afc) diXjdjXk — {Xj — Xk) djXidiXk 

The semihamiltonity (1.4) then is written as the condition of vanishing the quan- 
tities (4.7) 

(4.8) al,^ =0, for z ^ fc ^ j 
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Partial derivatives in (4.7) as in (1.4) are calculated with respect to the vari- 
ables u^, . . . ,u" for which the matrix A(u) is diagonal. The frame of eigenvectors 
X^, . . . , X" is chosen to be the coordinate frame for such variables. 

Using the Froelichcr-Nijenhuis bracket we construct the tensor K from the ma- 
trix A as follows 

(4.9) K = 3{{A,A},A2} = 3{N,A2} 

Tensor K defines the vector- valued 3-form the values of which for the frame vectors 
can be calculated according to (4.6). As a result from (4.9) we have 

(4.10) K(Xfe, X„ X,) = Kl^^u + i^^,,X, + Kl^y., 

(there is summation by i,j, fc here). We are interested only in one group of com- 
ponents of tensor K. Others can be obtained by cyclic transpositions of indices 
in if^j,. The values of components from this group in (4.10) are defined by the 
following formula 

^4 ^^^ Kl, " aL, -2(A, - A,) [(A, - A,) + (A, - A,)] 

{diXkdjXk — diXjdjXk — djXidiXk) 

From (4.11) we see that the difference K^^- — a^j ■ contains only the derivatives of the 
first order. Now taking the tensor N we construct another tensor M. Corresponding 
polylincar form is the following 

M(X, Y, Z) = N(X, AN(Y, Z)) + 
+ N(AX, N(Y, Z)) - N(N(X, Z), AY)-|- 
^ '^ ' +N(N(X,Y),AZ) -N(X,N(AY,Z))- 

-N(X,N(Y,AZ)) 

For the tensor M from (4.12) computed in the frame of eigenvectors of the operator 
A(u) we have 

(4.13) M(X,, X„ X,) = M,%-Xfe + M^,^-X, + M^.^.X, 

For the of components components of M in (4.13) we get 

Mfc% = -(A«-A,)(A.-A,)(A, -Afe) 



(4.14) 

{diXkdjXk — diXjdjXk — djXidiXk) 

The coefficients M|j and M^-- aren't of interest for us now. Comparing (4.11) 
with (4.14) we define another tensor Q 

Q(X, Y, Z) = K(AX, AY, Z) - KiA^Ji, Y, Z)- 
(4.15) - K(X, AY, AZ) + K(AX, Y, AZ) + 4M(AX, Y, Z)- 

- 2M(X, AY, Z) - 2M(X, Y, AZ) 
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For the components of this tensor in the frame of eigenvectors of operator A(u) we 
get the relationship 

(4.16) Q(Xfe, X„ X,) = Ql^Xk + OL,X. + QL,X, 

which is analogous to (4.10) and (4.13). The components in (4.16) which are of 
interest for us are expressed through (4.7). They have the form 

(4.17) Ql,^ ^ -{X, - Xk)iX, - Xk)ai, 

Now on a base of (4.17) we are able to construct the semihamiltonity tensor which 
is the main goal of the whole paper. It is defined by the following formula 



P(X, Y, Z) = AQ(X, AY, Z) + 

+ AQ(X, Y, AZ) - A2Q(X, Y, Z) - Q(X, AY, AZ) 



(4-18) , .^.^.. ... ,2 



It is easy to check that for arbitrary three vectors from the frame X^, . . . , X" one 
has the relationship 

(4.19) P(Xfc,X„X,) = (A, - XkfiXj - Xkfal^^Xk 

As a result we proved the following theorem (the invariant criterion of semi- 
hamiltonity) . 

Theorem 9. The diagonalizable operator of general position A with mutually dis- 
tinct eigenvalues is semi-Haniiltonian if and only if when associated tensor P from 
(4.18) is identically zero. 

The proof of this theorem follows directly from (4.8) and (4.9). It doesn't require 
any comments. Concluding all above considerations we give the formulae which 
enable us to calculate the components of the tensor of semihamiltonity P from the 
matrix A(u) 

n n 

This formula is derived from (4.18). Components of the tensor Q in the formula 

(4.20) are calculated on a base of (4.15) 

n n 

^kij — /_^ Z^y k^pqj^i + ^k^piq^j ^ ^q^k^pij~ 
P=l «=1 

n 

- i^|,,Af AJ) +J2{AAIm;,^ - 2M4^.Af - 2M,%A^) 
p=i 

Components of M in the above formula are found from (4.12) 

n n 
p=l q=l 
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Tensor K are computed through tensor of Nijenhuis with the use of the bracket of 
FroeUcher and Nijenhuis on a base of formula (4.9). Let's take B = A^. Then for 
the components of K we have 



^4 21) K^,,,^Y.iBp9kNf^-Bld,N., 



p=i 






By dots in (4.21) we denote 10 summand that can be obtained from 5 cxpHcit 
summands by cychc transposition of indices i, j and k. 

The above formulae huge enough for direct calculations. However modern com- 
puter systems for analytical calculations solve this problem for any particular equa- 
tions of hydrodynamical type in applications. 
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